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I. INTRODUCTION 



The notion of the entanglement in bosonic or fermionic systems has been the subject 
of many recent discussions and controversy since it rises some conceptual questions about 
the nature of the quantum correlations appearing in such systems[l-5|. The main difficulty 
appearing in the definition of a criterion apt to classify and quantify the amount of quan- 
tum correlations is due to the intrinsic indistinguishability of the particles supposed to be 
entangled, which results in the symmetrization or antisymmetrization of the wavefunction. 
This stimulated great interest in the scientific community and lead to various proposals of 
entanglement criteria for identical particles. In particular, a number of works present in the 

n nti n 

literature have been focused on the analysis of bipartite entanglement [l|, |3|, |4|, |6j. In the 
paper by Schliemann et al. [1], the amount of quantum correlations between two fermions 
is evaluated by considering the fermionic analogous of the Schmidt rank, namely the Slater 
rank, which is given by the number of Slater determinants needed to expand the entangled 
states. The Schliemman's criterion has been used to evaluate the entanglement in a number 
of bipartite systems of physical interest [7H9| , still it is discussed since it does not seem to 
behave correctly under one-site (local) and two-site (non local) transformations [2j. In the 
approach developed by Zanardi [3] , the entanglement is evaluated in terms of the quantum 
correlations between modes by mapping the Fock space of the modes themselves into qubit 
states. It has been shown that such a criterion overestimates the entanglement between the 
two parties of the quantum system since it conflicts with the local number particle superse- 
lection rule \&\. In order not to violate the latter, Wiseman and Vaccaro J4] have proposed an 
operational criterion for bipartite entanglement of identical particles. They define as entan- 
glement of particles the maximum amount of quantum correlations which can be extracted 
by means of local operation on the modes from the two parties of the system and then set 
in standard quantum registers (composed of distinguishable qubits) 0, 6]. Only recently a 
coherent picture on the entanglement of indistinguishable particles somehow including all 
the above criteria has been presented in the literature |5(. 

Though the entanglement of multipartite systems suggests more possibilities and phenom- 
ena with respect to the bipartite case, so far it has not been widely investigated. Its quan- 
tification would undoubtedly represents a key ingredient to understand a number of physical 
phenomena involving correlated many-particle systems such as quantum phase transitions, 
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quantum Hall effect, high temperature superconductivity 1CH12I]. Among the various phys- 
ical systems of identical particles, particular attention has been devoted to the analysis of 
multipartite entanglement in a non-interacting fermion gas. It has been shown that the 
multiparticle entanglement can be built only out of two-fermion entanglement [3]. The 
existence and the nature of the entanglement in systems of noninteracting three fermions 
have also been analyzed by means of different parameterized entanglement witnesses [l^-[l6] 
that enable to detect two different classes of tripartite entangled states,namely Greenberger- 
Horne-Zeilinger (GHZ) and W states. By using the Zanardi approach [3J, that is by mapping 
the Fock space of a systems of fermions into the isomorphic mode state, a geometric mea- 
sure of entanglement for N fermions with spin ~ has also been proposed [171 ] . A general 
entanglement criterion applicable to both bosonic and fermionic systems has recently been 
introduced though it is defined only with reference to the measurement setup 

In this work, we discuss an entanglement criterion suitable both for three-fermion and for 
three-boson mixed states. We follow the basic concepts of the method proposed by Wiseman 
and Vaccaro [4| for the analysis of the bipartite entanglement in systems of identical particles. 
Our criterion does not violate the particle local number superselection rule. In fact here 
the entanglement between three subsets, say A, B, C is given by the average, over the 
local particle number, of the amount of quantum correlations existing among the standard 
quantum registers of three subsets each with a definite local particle number. Given the use of 
concept of standard quantum registers, the tripartite negativity (TPN), recently introduced 
to quantify the tripartite entanglement of non identical subsystems Q], is here taken as 
a measure of the entanglement exhibited by three-boson or three-fermion mixed states. In 
order to get a better understanding of the criterion, first we compare and contrast it to other 
approaches, evaluating the quantum correlations in many-particle systems. Then we apply 
it to investigate the entanglement dynamics in a simple model mimicking the continuous- 
time quantum walks (QWs) 20|] of three identical particles (bosons or fermions) in an one- 
dimensional structure. Specifically, QWs of indistinguishable elements constitute the ideal 
systems to test and validate our theoretical approach since in them the exchange terms 
lead to the formation of correlations between the particles even in absence of interaction, as 
shown by the experimental investigations on two-photon propagation in lattices of coupled 
waveguides |2l|, [22 ] . 

The paper is organized as follows: In Sec.[TT]we define the tripartite entanglement criterion 
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in boson and fermion systems with reference to other entanglement measures in systems of 
identical particles. The main properties of the TPN as a measure of the amount of quantum 
correlations in a six-mode system and its comparison with the geometric measure introduced 
in Ref. are discussed in Sec. IIHI In Sec. IIVI we evaluate the time evolution of the TPN 
in QWs of three particles propagating in a lattice for two different cases: particles following 
the Bose-Einstein statistics and particles following the Fermi-Dirac statistics. Finally, we 
present our conclusions in Sec. [V] 

II. TRIPARTITE ENTANGLEMENT CRITERION 

In this section, we introduce an entanglement criterion for boson and fermion systems 
relying on the concept of entanglement of particles [4j . 

In order to represent the mixed quantum state of N identical particles, we adopt the 
occupation-number representation. Thus, in a system with L modes an arbitrary mixed 



state of iV-bosons (fermions) can be written as 231 ]: 



P= E /(W,K})i{n})({n'}| (1) 
{«}>'} 

with the integers (n-) of the set {n} = m, . . . , . . . , ni ({n 1 } = n^, . . . , . . . , n' L ) 
satisfying n\ . . . + Hi . . . + n L = N (n' x . . . + v! i . . . + n' L = N) . The ket \{n}) is the state of 
the Fock space with rij's particles in the single-particle modes z's and f({n}, {^'}) are the 
coefficients in the superposition. While for bosons n^s range from to N, for fermions the 
occupation numbers are restricted to and 1 by the Pauli exclusion principle. 

A formal equivalence between the Fock space and the space obtained as tensor product 
of L-dimensional single-particles subspaces can be established [3J. In such a procedure, the 
occupation number of each mode represents a distinct state of the mode itself. Zanardi, 
Lari et al. 0, [l^] defined as entanglement of a system of identical particles the so-called 
entanglement of modes (Em), namely the amount of quantum correlations appearing among 
the occupation numbers of the modes. 



As argued by some authors 



24 



251 ] . an application of the E M , as it is, to the quantum 
states of many-body systems could lead to misleading results since such a measure quantifies 
the entanglement between modes and not between the particles. In order to clarify this 
point, let us consider a single particle in an equal superposition of three modes, so that in 
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the occupation- number representation its state \x) reads 



\ X ) = -j= (|001) + |010) + |100». 



(2) 



By using the geometric measure of the entanglement relying on the isomorphism between 
the Fock space and the mode space 17| . we can evaluate the tripartite entanglement tQ 
between the modes as: 



T 



sep; 



where 



T 



a Kxki ® o-j ® ^fclx)! 2 , and ||r| 



sep 



1 



(3) 



(4) 



with {<7j} the set of three generators of the SU(2) group, namely the Pauli operators, each 
acting on a single-mode space. After straightforward calculations, one finds that the tri- 
partite entanglement of \x) is different from 0, specifically ec= a/33/3 — 1. Such a result 
would seem to indicate a single-particle non-locality stemming from the entanglement with 



the vacuum 



24 



251 ] . However, the above entanglement is not apparent when the particle 



wavefunction x( x ) is considered in the space configuration 



(5) 



with ipi indicating the single-particle states. 

To explain the emerging paradox, the appearance of non-locality has been related to 
multi-particle effects 25j. Anyway, this behavior arises some questions on the use of the 
criteria relying on the E u , since they do not always capture the true entanglement between 
the parties of the system. As argued by Wiseman and Vaccaro [4j, in particular such criteria 
fail to take into account the local particle-number (LPN) superselection rule. In fact each 
party of the system must be able to perform arbitrary local operations on its modes in order 
to fully use the entanglement between the modes. Unless each subsystem possesses a definite 
number of particles, such local operations violate the superselection rule for the LPN and 
are not possible in practice. Thus in general, the entanglement measures relying on Em 
overestimate the available entanglement. 

A tripartite entanglement criterion obeying the local particle-number superselection rule 
can be obtained by extending the operational definition of bipartite entanglement of identical 
particles introduced by Wiseman and Vaccaro j^, 6]. Being Alice, Bob, and Charlie the 
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three parties of a quantum system in the mixed state p each accessing to a given set of 
modes, here we assume that any subsystem possess a standard quantum register, namely 
a set of distinguishable qubits in addition to the indistinguishable particles described by 
p. We define as tripartite entanglement er the maximum amount of quantum correlations 
that Alice, Bob, and Charlie can produce between their standard quantum registers by 
means of local operations. As a consequence of the LPN superselection rule, this tripartite 
entanglement in place of quantum correlations between the modes that Alice, Bob, and 
Charlie have access to, is given by 

n A +n B +n c =N 

e T = 

Pn A ,n B ,n c e ABc{Pn A ,n B ,n c ), (6) 

n A ,n B ,nc 

where p nA ,n B ,n c = ^-n A ,n B ,n c P^-n A ,n B ,n c is obtained from p by means of the projectors 
Tl nA ,n B ,n c onto fixed LPN's states (n^ for Alice, % for Bob, and n c for Charlie). P nA ,n B ,n c = 
Tr (pn A ,n B ,nc) * s the probability of finding riA,n B , and n c as a result of a measurement of 
the local particle number by Alice, Bob, and Charlie, respectively, while €abc is an entan- 
glement standard measure which quantifies the degree of tripartite entanglement among the 
three sets of modes each controlled by a party of the system. It is worth noting that here is 
used a standard measure of the tripartite entanglement, that is a measure of the tripartite 
entanglement of non identical subsystems, since the standard quantum registers of Alice, 
Bob and Charlie consist of distinguishable qubits. Specifically, since some good measures 
tAB of bipartite entanglement can successfully be extended to multi-partite systems by 
considering bipartite partitions of them 2g|, we take as standard tripartite entanglement 
measure the geometric mean of the entanglement measures of the bipartitions of the system: 

^ABC — \J ^A-BC e B~AC^C-AB (7) 



in agreement with what proposed in the literature 18|, Q, [2^ . Thus, expression ([6]) now 
reads 

n A +n B +nc=N 



e T — Pn A ,n B ,n c y e A-Bc(Pn A ,n B ,n c )^B-Ac(Pn A ,n B ,n c ) e C-AB(Pn A ,n B ,n c )- (&) 

n A ,n B ,n c 

The tripartite entanglement €t does not violate the LPN superselection rule. In fact 
when Alice, Bob, and Charlie measure their local particle numbers ha,ub, and nc, the 
mixed state p collapses, with a given probability Pn A ,n B ,n c -> m Pn A ,n B ,n c - Local operations 
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can be performed on the latter without any restrictions in order to transfer its entanglement 
e(Pn A ,n B ,nc) to the standard quantum registers of each subsystem. Since et is the weighted 



sum of the terms e A - BC (p nA ,n B , nc ) e B- ac (pn A ,n B ,n c ) e c-AB {Pn A ,n B ,n c )i the entanglement is 
not affected, on the average, by measurements on the local particle number. 

Unlike Eq, the tripartite entanglement of the single-photon state written in Eq. (j2J) ac- 
cording to our criterion is 0. By considering the partition of the system in three subsystems 
each controlling a single-mode state, the state \x) is given by the linear superposition of three 
quantum states having different LPN's and therefore the contribution to et from each terms 
is zero. This result again supports the interpretation of et as a measure of the tripartite 
entanglement not among the modes but among the particles of the systems, in agreement 



with the analogous investigations performed for the bipartite case 



III. TPN IN A SIX-MODE SYSTEM 

In order to obtain et according to the criterion given in Eq. (JBJ), the evaluation of the 
bipartite standard entanglement ei~jK(Pn A ,n B ,n c ) f° r & U bipartitions I — JK of each possible 
set of LPN's is needed (with I— A, B, C and JK=BC, AC, AB, respectively). Unfortunately, 
the latter is a very challenging task in high- dimensional systems for mixed state of multi- 
particle systems. 

In this section, we intend to give a non ambiguous measure of the tripartite entanglement 
et which can be used in a practical way in a simple system and then compare it with the 



geometric measure [l7j. To this aim, we restrict our analysis to the case of three bosons 
(fermions) each in a six-mode single-particle space Hq. The Fock space Hq(3) of the full 
system is the totally symmetric (antisymmetric) subspace of the tensor product hf 3 . By 
partioning the system in three subsets, each acceding two modes, we note that only a few 
quantum states belonging to the Fock space Hq(3) give a non vanishing contribution to 
the tripartite entanglement according to the criterion given in Eq. (jSJ). Let us examine 
the case of the set of mixed states F( nA=0 ) := span{p( 0)riBiric .))} (riB,n c = 0,1,2,3 with 
n B + n c — 3) where the local number of particles possessed by Alice, is zero. Any state of 
F(n A= o) can be written as |00)(00| A <g> psc an d therefore it is biseparable, that is, it can be 
factorized in a term describing Alice with no particle and in a term describing the Bob and 
Charlie subsets with three particles. The bipartite entanglement relative to the bipartition 
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A — BC is zero and, in turns, is vanishing according to Eq. (jSJ) . This implies that the 
entanglement due to sets of quantum states with a LPN equal to zero is vanishing. On the 
other hand, the tripartite amount of quantum correlations can be different from zero only 
when each party has one particle, that is only for the quantum states belonging to the set 
span{p(i i i i i)}. Since each partition has two modes, the subspace spanned 
by the vectors of F( na=liri6=1) „ c=1 ) turns out to be isomorphous to the three-qubit Hilbert 
space C 2 <S> C 2 <S> C 2 . Thus the entanglement of three identical particles (both bosons and 
fermions) in this case can be evaluated in terms of the quantum correlations between three 
qubits. 

Different approaches have been developed to characterize the tripartite entanglement 
in system of distinguishable qubits. The three-tangle or residual entanglement related to 



28 



even if its use is questioned 



squares of the bipartite concurrences has been proposed 
since it does not properly quantifies the three-party entanglement for W states |29j. Alterna- 
tively, as stated in the previous section, standard tripartite entanglement can be estimated 
by means of some good measures of the amount of quantum correlations in bipartite systems 
whose Hilbert space is C 2 ® C 4 . It should be noticed that in this case both von Neumann 
entropy and concurrence presents weaknesses [3]. The former is an appropriate measure 
only for pure states, the latter, even if it is well defined for mixture states of two qubits, is 
applicable to higher dimensions only for pure states. A valid measure of tripartite entan- 
glement for non pure-state can be obtained by using the negativity. From the Eq. (EJ), the 
tripartite negativity (TPN) Mabc can be defined as 



Mabc = ZJNa-bcNb-acNc-ab , (9) 

where the bipartite negativities are given by Nj_ JK = J2i\lii.P Tl )\ ~ 1 with 7i(p T/ ) the 
eigenvalues of p Tl which is the partial transpose related to the subsystem / of the total 
density matrix. 

By inserting the TPN of Eq. (Q in Eq. (jSJ), the tripartite entanglement of a mixed state 
p of three bosons (fermions) in a six-mode system, can be easily evaluated from 



£t = Pi^y J^a-bc{piaaWb-ac{pi,i,iWc-ab{pi,ia), (10) 

once is known the projection of p on the subspace where the local number of each party is 
fixed to 1, namely ^l^n^ipll^i. 
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In order to illustrate some features of the entanglement measure introduced in Eq. (flOj) 
in comparison with the geometric measure 6q, here we consider a six- mode three-fermion 
system in the state 

sin a sin at 

|$) = cosacos/3|010101) + cos a sin /3| 101010) + — ^|111000) + — H000111) (11) 

v2 V2 

where a, j3 are phases ranging from to it. If we consider the following partition of the 
system A={1,2}, £?={3,4}, and C={5,6}, we note that |$) is a linear superposition of 
states, each with a given LPN for any set of modes. When a=0,n and P=(tt/A), (37r/4), 
only the first two terms are non vanishing and a GHZ state with each party of the system 
containing a particle, namely, nA=riB=nc=l, is obtained. On the other hand for a=(7r/2), 
a specific LPN cannot be ascribed to the subsystems. Our goal is to evaluate the amount 
of quantum correlations of |$) according to et and €q- The latter can be calculated again 



from the Eq. ([3]) , where now 17| : 



15 



J2 8|($|A 4 ® A^Afcl^)! 2 , and ||r|| sep = 6^6, (12) 



with {Aj} indicating the set of fifteen generators of the SU(4) group [30J acting on the two 
sites of one subset. 

Fig. [1] displays and as a function of a and 0. The two entanglement measures 
show a different behavior and in agreement with the theoretical predictions of the previ- 
ous section the geometric measure is always greater than the TPN. The latter is vanishing 
for any value of /3 when a = 7r/2; in this case, |<&) becomes the linear superposition of two 
terms, each corresponding to different LPN's thus yielding zero entanglement. At a=0,7r and 
/3=(7r/4), (37r/4), €t exhibits four peaks indicating the maximum amount of quantum corre- 
lation in the GHZ states l/y/2 ( ± |010101> ± |101010)^ describing one particle in each par- 
tition. On the other hand, the geometric measure is zero only for those values of the phases 
reducing the linear superposition of the four states in Eq. (II ip to a single term. For a = tt/2, 
€q shows a local minimum but unlike e-r it is not vanishing thus indicating a degree of non- 
separability different from zero. Its maximum values are found when /3=(7r/4), (37r/4) and 
a=(7r/4), (37r/4), namely when all the moduli of the coefficients of the linear superposition of 
Eq. dTJ]) become equal and |$) reduces to 1/2 ( |010101) ± |101010) ± |111000) ± |000111) j . 



The disagreement between the qualitative behavior exhibited by the tripartite negativity 
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FIG. 1. (Color online) The tripartite entanglement ec (left panel) and ex (right panel) as a function 
of the phases a, and f3 ranging from to tt. Specifically, the quantum correlations exhibited by 
the three-fermion state |<3?) of Eq. (|lip among the three parties, t4={1, 2}, B={3, 4}, and C={5, 6} 
are here considered. 



and the geometric measure for this specific case is representative of the discrepancy between 
the two entanglement notions. In fact, €q quantifies the amount of quantum correlations 
between the modes regardless the number of particles in any subset, while €t measures the 
entanglement of the particles in the modes. 
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IV. ENTANGLEMENT DYNAMICS IN THREE-BOSON- AND 

THREE-FERMION-QWS 

In the last years, theoretical and experimental investigations have shown how the quan- 
tum statistics, due to exchange symmetry, may play a key role in the appearance of quan- 
tum correlations between non-interacting identical particles in various systems, ranging from 



the propagation of photons 3lJ to electron transport in integer quantum Hall effect 32]. 
Among the different physical phenomena examined, the continuous-time quantum walks 
(QWs) have received great interest since they represent an ideal laboratory to observe 
many-particl e q uantum mechanical behavior and to implement future quantum technolo- 



gies 



21 



22 



33 



341 ] . Specifically, the emergence of non classical correlations in QWs of two 



photons in an array of coupled waveguides has experimentally been observed [2l|, |22J. By 
using the criterion based on TPN described in the previous section, here we investigate the 
entanglement dynamics of the continuous-time QWs of three bosons and of three fermions 
in an array of six sites. 

In order to study the three-particle entanglement dynamics, we follow the formalism 



developed in Ref. [33| for the QWs of non-classical light in an array of coupled waveguides. 
We consider a one-dimensional tight binding model, within the approximation of nearest- 
neighbor interaction, described by the Hamiltonian H: 

6 5 

U = G 4* + T J>tc w + 4 +lCi ), (13) 

i=l i=l 

where cf (q) is the creation (annihilation) operator for a particle in the i-th site which 
obeys the commutation relation [q, c^]±=6ij, where the upper sign holds for bosons and the 
lower one for fermions [35]. G is the on-site energy and T is the tunneling rate at which 
the particles are transferred to the neighboring sites. Here we are considering a lattice with 
reflecting boundary conditions, so that the transition amplitude from the node 1(6) to the 
node 6(1) is zero. As it will be shown in the following, this implies that the time evolution 
of the quantum system is not periodic unlike the QWs of particles in lattices with periodic 
boundary conditions 0, Q • 

By diagonalizing the Hamiltonian ffl3|) . in the Heisenberg picture the annihilation oper- 
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ator c r at time t can be obtained as 33[: 



c r(^) = y^c 5 (o)c rs , 

8=1 

C rs = ^ exp ^ — — j ^ exp ^ — cos J j sm ^— j sm ^— j , (14) 



fe=i 

^6 



where | 5^fc=i srn (^y 1 ) sin ( § ^ L )—^rs- Since any input Fock state can be expressed by means 
of the creation operators cj an d the vacuum state |0, 0, 0, 0, 0, 0), its time evolution during 
the QWs can be calculated by using the Eq. f fT4|) . Even if the effect of different initial 
occupation numbers on the entanglement dynamics could be considered, here, for simplicity, 
we only examine the case of three particles (bosons or fermions) initially coupled to three 
neighboring sites, i.e. the state 

|M>(t = 0)> = |1, 1, 1, 0,0, 0>. (15) 

By means of the time evolution of the quantum state |\l/(t)), we can quantify the tripartite 
entanglement dynamics of the system by using the TPN defined in Eq. f JTQj) . To this purpose, 
a numerical approach has been used which allows one to evaluate at any t from |^(t)) the 
density matrix pi,i,i where the LPN of each subsets is fixed to 1 and then to diagonalize 
its partial transposes for the evaluation of the bipartite negativities. Here we consider two 
different partitions of the system, namely A={1,2}, £>={3,4}, C={5,6} and yl'={l,4}, 
£?'={2,5},C"={3,6}. 

In Fig. [2] we report the entanglement of the three-boson and the three- fermion state as a 
function of the "time" r = tT/h for the two above partitions. We note that in all the cases 
the amount of quantum correlations in the tripartite system under investigation oscillates 
with time though such oscillations are not periodic. In fact the interference of the three- 
particle wavefunction amplitudes, due to the backscattering at reflecting boundaries, results 
into a non-periodic time evolution of the system. Unlike the case of single-particle dynamics 



in QWs into periodical structures 



36. 



371 ]. a time revival, namely the time interval needed 
to reconstruct the whole initial wavefunction, cannot be found here. 

For both of the partitions, we note that the tripartite entanglement for the fermions 
is always greater than the one exhibited by the bosonic systems. This behavior can be 
explained by taking into account the quantum statistics of the particles. The Bose-Einstein 
particles exhibit bunching thus making finite the probability to find more particles in the 
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A={1,2},B={3,4},C={5,6} A'={ 1,4}, B'={2,5},C'={3,6} 




X X 

FIG. 2. (Color online) Entanglement as a function of the "time" r for the fermionic (blue line) 
and bosonic (green line) systems when two different partitions are considered: A={1, 2}, B={3, 4}, 
C={5,6} (left panel) and A'={1,4}, £'={2,5}, C'={3,6} (right panel). 

same site or in the same subset during the time evolution of the system. On the other hand, 
due to the exclusion Pauli principle, no more than one fermion can occupy the same mode. 
This implies that the term P\,i,i describing at any time the probability to find one particle 
in each subset and appearing in the definition of TPN of Eq. ( ITU]) , increases moving from the 
boson- to the fermion-system, even if both of them are initially described by the same state 
\*$!(t = 0)). As a consequence, the tripartite entanglement among bosons results smaller 
than the one among fermions. 

The TPN dynamics of the two partitions examined presents some differences. While 
the tripartite bosonic entanglement shows the same qualitative behavior, that is oscilla- 
tions within the interval ranging from to 0.2, the amount of quantum correlations among 
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FIG. 3. (Color online) Single and two-particle density probability of the bosonic and fermionic 
systems initially described by the state |^) of Eq. (|15p for r=8.7. (a) The single particle density 
distribution p r which is identical for fermions and bosons, (b) The correlation matrix T rs for bosons 
computed from Eq. f)16[) . Bosons tend to localize in the the same or different sites of the first half 
of the lattice, (c) The correlation matrix T rs for fermions computed from Eq. f)lT[) . Fermions 
are found with high probability in the neighboring sites of the first half of the lattice, (d) The 
interparticle distance probability for bosons. They tend to appear at same site or at neighboring 
sites, (e) The interparticle distance probability for fermions. They are more likely to be separated 
by a small number (1-2) of sites. 



fermions in the partition where subsets are composed of non adjacent modes results, on the 
average, larger than the one found for subsystems controlling neighboring sites. In order to 
deeper investigate this behavior, we focus on some single- and two-particle features of the 
quantum system at "time" r = 8.7 where the fermionic TPN between A', B', C shows a 
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peak. In Fig. [3] the single particle density p r {t) = (cjc r ) and the two-particle correlation 
F rs (t) = {clc\c s c r ) are reported. Here we do not consider the three-particle correlation func- 
tion, since the two-particle one is sufficient to understand the features of the entanglement, 
as will be shown in the following. p T {t) does not permit to analyze the quantum properties 
of the system under investigation since it shows the typical features of the classical QWs of 
two uncorrelated particles. As a matter of fact, the single-particle density whose expression 
in terms of the coefficients C rs reads p r (t) = X/s=i I^Vsl 2 ^ (with n s the number of particles 
in the site s at t=0) does not depend upon quantum statistics and turns out to be identical 
for bosons and fermions (see panel (a) of Fig. [3]). 

On the other hand, in analogy with the Glauber photodetection theory [38[ T ra (t) indi- 
cates the probability to detect at time t one particle at site r and the other one at site s and 
displays the two-particle quantum correlations between the locations r and s depending on 
the quantum statistics of the particles. For bosons it is given by 

6 p-l 6 
Trs^) = ^ ^ \C rp C S q 4~ C r gC S p \ 2 n p n q + \C rp \ 2 \C sp \ 2 n p (n p - 1), (16) 

p=l q=l p=l 

while for fermions it reduces to 

6 p-l 

Trs(^) ^ ^ ^ \C rp C S q CfqC ' S p\ f^plT'q- (1^) 
p=l q=l 

As shown by the correlation matrix in the panel (b) and (c) of Fig. [31 bosonic bunching 
results into a very large probability to find two particles in the same sites in the first half of 
the lattice. Instead if the particles follow the Fermi-Dirac statistics, they tend to separate 
to neighboring sites of the first half of the lattice (see panel (c) of Fig. [3]). Regardless the 
quantum statistics, this implies that for the partition {A, B,C} the probability that more 
particles occupy the same subset is very high and, in turns, this means a low value of the 
entanglement as depicted in the left panel of the Fig. [2j On the contrary, for the partition 
where the subsystems control nonadjacent modes, the LPN of A',B', and C is more likely 
to be equal to 1 and therefore the TPN can reach high values (see the right panel of the 
Fig. E]). 

Also the interparticle distance probability [34] g(A) = ^2 q ^ q , q +A, depicted in the panels 
(d) and (e) of Fig. |3] supports our interpretation. Apart from the behavior in A=0, related 
to the role played quantum statistics in determining the probability to find two particles in 
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the same site, g is high for low values of the intersite distance both for bosonic and fermionic 
system thus confirming the localization of particles in neighboring sites. 

A more exhaustive discussion of the entanglement dynamics in our model undoubtedly 
requires further analyses also on three-particle quantum properties of the system at any time 
but this goes beyond the scopes of this work. 



V. CONCLUSIONS 

In summary, we have presented a criterion for the tripartite entanglement of mixed states 
of indistinguishable particles (both bosons and fermions). The entanglement measure pro- 
posed deviates from the other approaches recently advanced [17] . They rely on the entan- 
glement of modes Q| but their physical significance is not always very clear as shown here by 
some misleading results. As a matter of fact, such measures can give values different from 
zero even when applied to single-particle states which describe a particle in equal superpo- 
sition of modes. In this paper, we developed a different criterion based on the notion of 
bipartite entanglement of particles as introduced by Wiseman and Vaccaro [4j . It provides 
a valid guideline to quantify the amount of quantum correlations in tripartite systems on 
which the allowed physical operations, such as preparation, manipulation, measurement, do 
not change the local number of particles contained in each subset. 

Specifically, we have shown that a good entanglement measure satisfying the requirements 
of the theoretical criterion introduced in this pap er can be obtained from the TPN recently 
used in the context of non identical particles [lfl, Q . Indeed, it constitutes a very practical 
tool to evaluate the amount of quantum correlations appearing in three-boson or -fermion 
mixture states of low- dimensional systems, that can be extracted and then placed in standard 
quantum registers without any violation of the local particle number superselection rule. As 
an application, we have considered a simple toy-model of three fermions in a system with six 
modes and proved the realibility of the TPN in comparison with the geometric multipartite 
entanglement measure introduced in Ref. [171. In agreement with theoretical predictions, 
we found that TPN is greater than the geometric measure. The qualitative behavior of the 
two criteria as a function of the parameters of the model is quite different thus indicating 
their conceptual discrepancy. The geometric measure gives an estimation of the amount of 
quantum correlations among the occupation numbers of the sites regardless the particles. 
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On the other hand the TPN, even if it depends upon the local particles number in the 
subsystems, quantifies the entanglement among the particles themselves. 

Furthermore, we have used the TPN to quantify the time evolution of the quantum 
correlations for the continuous-time QWs of three non interacting particles (both bosons 
and fermions) in a six-mode lattice, for two different partitions of the system. As expected, 
we find that the entanglement dynamics is not periodic due to the presence of reflecting 
boundary conditions in the lattice sites. Furthermore, it strongly depends upon the quantum 
statistics of the particle involved in the process: the bosonic bunching leads to a large 
probability to find more bosons in the same site and therefore in the same subset of the 
system. This implies a local particle number equal to zero for another of the two parties and 
therefore the separability of the three-particle state. On the other hand, for fermions the 
Pauli exclusion principle prevents occupation number of sites greater than one, thus making 
more likely to find correlated states with one particle in each subset. 

Although the criterion introduced in this paper has been mainly applied to the investiga- 
tion of the tripartite entanglement in some simple models, such as continuous-time QWs, it 
is completely general, and its application turns out to be helpful to investigate the building 
up of quantum correlations in both bosonic and fermionic systems. In particular, it would be 
of interest to extend the TPN to those physical system with large dimensionality where new 
and interesting results could be expected due to higher complexity of the systems themselves 
and to consequent richer structure of the quantum correlations. 
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